Nonequilibrium Transport through a Kondo Dot: 
Decoherence Effects. 
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We investigate the effects of voltage induced spin-relaxation in a quantum dot in the Kondo 
regime. Using nonequilibrium perturbation theory, we determine the joint effect of self-energy and 
vertex corrections to the conduction electron T-matrix in the limit of transport voltage much larger 
than temperature. The logarithmic divergences, developing near the different chemical potentials of 
the leads, are found to be cut off by spin-relaxation rates, implying that the nonequilibrium Kondo- 
problem remains at weak coupling as long as voltage is much larger than the Kondo temperature. 
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Electron transport through quantum dots or point con- 
tacts possessing a degenerate ground state (e.g. a spin) 
is strongly influenced by the Kondo effect- , provided the 
dot is in the Coulomb blockade regime. In the linear re- 
sponse regime, the Kondo resonance formed at the dot at 
sufficiently low temperature, i.e. at or below the Kondo 
temperature Tk, allows for resonant tunneling, thus re- 
moving the Coulomb blockade and leading to conduc- 
tances near the unitarity limit. This has been observed 
in various experiments on quantum dot devices^. 

The Kondo resonance is quenched by either large tem- 
perature, T ^> Tk, large magnetic field, B » Tr or a 
large bias voltage, V ^> Tk- However, the mechanism of 
how and why the Kondo effect is suppressed is qualita- 
tively different in the three cases. The Kondo effect arises 
from resonant spin-flip scattering at the Fermi energy. 
Temperature destroys the resonance mainly by smearing 
out the Fermi surface, whereas a magnetic field lifts the 
degeneracy of the levels on the dot and thereby prohibits 
resonant scattering. The effect of a bias voltage, V, is 
more subtle. It induces a splitting of the Fermi energies 
of the left, and the right lead. However, this splitting af- 
fects directly only resonant electron scattering from the 
left to the right lead, but not any scattering which begins 
and ends on the same lead. Yet these remaining resonant 
processes are suppressed by a different effect: the voltage 
induces a current which leads to noise and therefore to 
decoherence of resonant spin- flips. It is the goal of this 
paper to study those decoherence effects in detail. 

In perturbation theory, the signature of Kondo physics 
are logarithmic divergences arising from (principle value) 
integrals of the type 
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where f(u>) is the Fermi function, D a high energy cutoff 
(i.e. bandwidth) and Em some infrared cutoff. There are 
three rather different ways to cut off the logarithm, and 
to destroy the Kondo effect, corresponding to the three 
mechanisms discussed above. First, temperature broad- 
ens f(u>) leading to Em ~ T. Second, a magnetic field B 



shifts the pole with respect to the Fermi-energy, replacing 

h by ZJ^B' an d m this case ^IR ~ B. The third way to 
quench the logarithm is to introduce a finite decoherence 
rate r s , replacing £ by , implying E m ~ T s . 

The relaxation rate T s — T S (V, B, T) and the associated 
decoherence effects also exist in equilibrium. In the limit 
of vanishing bias voltage and magnetic field, the scale 
T s tends to a temperature dependent (Korringa) rate^, 
r s (0,0,T) < T, which vanishes as T — > 0, allowing for 
the quantum coherent Kondo state to be formed. In the 
case of a finite magnetic field and zero temperature, a 
B- and spin-dependent rate^, T Sirr (0,B,0) remains finite 
for the excited state a = j. In dynamic quantities it pro- 
hibits singular behavior at u> ~ B but it is not important 
for static quantities, where B eliminates all relevant sin- 
gularities. In the case of a finite bias voltage V, however, 
the finite rate r s (V, 0, 0) is instrumental to cut off singu- 
larities even in static quantities for T, B — > 0. The Kondo 
effect develops only to a certain extent, depending on the 
ratio V/T K . 

Not only for a quantitative description of experiments 
in the regime V 3> Tk, but even for a crude qualitative 
understanding of Kondo physics out of equilibrium, it is 
necessary to identify the correct relaxation rate T s . The 
question, how logarithmic contributions are cut off, is es- 
sential to derive the correct perturbative renormalization 
group description^ and to identify regimes where novel 
strong-coupling physics is induced out of equilibrium. 

The importance of the broadening of the Zccman levels 
was pointed out three decades ago by Wolf and Loseei 
in the context of the Kondoesquc tunneling anomaly ob- 
served in various tunnel junctions. Incorporating a Kor- 
ringa like, T and B dependent, spin-relaxation rate into 
Appelbaum's perturbative formula for the conductance^ 
was found to improve the agreement with experiments 
considerably (cf. e.g. Refs. |JJ and UXt) . Later, in the 
context of quantum dots, Meir et alm^ pointed out that, 
even at T = B = 0, the finite bias- voltage induces a 
broadening of the Zeeman levels. In their self-consistent 
treatment of the Anderson model, using the non-crossing 
approximation (NCA), this nonequilibrium broadening 
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was shown to suppress the Kondo peaks in the local den- 
sity of states, located at the two different Fermi levels. 
In Ref. ^3 we showed that this NCA relaxation rate is 
sufficiently large to prohibit the flow towards strong cou- 
pling for V ^S> Tk- In a perturbative study of the effects 
of an ac-bias, Kaminski et al. 13 argued that an irradia- 
tion induced broadening serves to cut off the logarithmic 
divergence of the conductance as T and V tend to zero. 
Treatments of the Kondo modeli^ and related problems^ 
at large voltages, which neglect the influence of deco- 
herence, find strong coupling effects even for V ^> T K . 
Coleman et ali^ recently argued that this is the case be- 
cause T s remains sufficiently small due to a (supposed) 
cancellation of vertex and self-energy corrections. 

To our knowledge, even to lowest order in perturba- 
tion theory, a systematic calculation of the nonequilib- 
rium decoherence rate is still lacking. It is the objective 
of this paper to provide such a calculation. This is a 
delicate matter since self-energy, and vertex corrections 
may indeed cancel partially, and an infinite resummation 
of perturbation theory is required. RecentlyiSiii, it was 
demonstrated that the Majorana fermion representation 
for the local spin- 1/2 circumvents this complication when 
calculating spin-spin correlation functions. In this repre- 
sentation, such correlators take the form of one-particle, 
rather than two-particle, fermionic correlation functions, 
and consequently only self-energy corrections have to be 
considered. Whether this representation will prove to 
be equally efficient for calculating other observables like 
the conduction electron T-matrix or the conductance re- 
mains to be seen. 

Based on the conjecture that no unexpected cancel- 
lations occur, we have recently developed a perturba- 
tive renormalization group description^ of the Kondo ef- 
fect at large voltages. In this approach, it was essen- 
tial to include the effects of T s . For usual quantum 
dots, the Kondo effect is sufficiently suppressed by TjtiH, 
such that renormalized perturbation theory remains ap- 
plicable at all temperatures, provided ln(V /Tk) 1. 
We argued that F,, as a physically observable quantity, 
should be identified with the transverse spin relaxation 
rate T2 = I/T2, measuring the coherence property of 
the local spin (More precisely, slightly different rates en- 
ter into various physical quantities, but to leading order 
in 1/]xi[V/Tk] one can use T s « ^2)- In this paper we 
show that within perturbation theory this is indeed the 
case, thus confirming our initial conjecture. Note that 
in more complex situations, for example in the case of 
coupled quantum dots, T s can be sufficiently small^ so 
that novel (strong coupling) physics can be induced for 
large voltages. 

In a preceding paper—, henceforth referred to as I, 
we calculated perturbatively the local magnetization and 
the differential conductance of a Kondo dot, including all 
leading logarithmic corrections in the presence of finite 
V and B. As effects of T s are not included to this order, 
some logarithms were not cut off by V but appeared to 
diverge with \n(D/T) or ln(D/\V — B\). A systematic 



calculation of the cut-off T s requires a consistent resum- 
mation of self-energy and vertex corrections. As will be- 
come clear in the following, this is a formidable task, and 
we have therefore concentrated on the quantity which 
appears to be most tractable: the conduction electron 
T-matrix as a function of frequency, in zero magnetic 
field. 

In Sec. [3 we introduce the model and some conventions 
used for the Keldysh perturbation theory. A combina- 
tion of self-energy corrections from Sec. Ill Al and vertex- 
corrections calculated in Sec. Ill Bl determines the spin- 
relaxation rate fSec. Ill 01) . In Sec. IIIII we show how this 
decoherence rate cuts off logarithmic corrections in the 
T-matrix. In Sec. HVI we consider the case of anisotropic 
exchange couplings and determine the exact combination 
of transverse and longitudinal spin- relaxation rates which 
enters the logarithms in the T-matrix. Appendices lAl and 
IbI contain details pertaining to Sections |ll Bl and lllll Ap- 
pendix El investigates how power-law singularities of the 
strongly anisotropic Kondo model are modified out of 
equilibrium by mapping it to the nonequilibrium X-ray 
edge problem for vanishing spin- flip coupling. 

I. MODEL AND METHOD 

We model the quantum dot by its local spin S (S = 5), 
coupled by the exchange interaction J aa ' (a, a' = L,R) 
to the conduction electrons in the left (L) and right (R) 
leads 

H = ^2 ( e k - Ma)c^ k(T c QkCT - g^ B BS z 

a,k,cr 

+ ^ ] Ja'a S ■ — C^, k , cr /T (T ' cr C Q i <;(T , (2) 
a, a' ,k,k' ,0.0' 

where Jlr describes a co-tunneling process transfer- 
ring an electron from the right to the left lead. Here 
Hl.r = ±eV/2 are the chemical potentials of respectively 
the left and right leads, r is the vector of Pauli matrices, 
g^sB the Zeeman splitting of the local spin levels in a 
magnetic field B, and c^ k(7 creates an electron in lead 
a with momentum k and spin a. We will use dimen- 
sionless coupling constants g aa > = N(Q)J aa ', with N(0) 
the density of states per spin for the conduction electrons 
(assumed flat on the scale eV, g^sB). We define for later 
use g d = {g LL + g RR )/2, Ag 2 = g\ L + g RR + 2g 2 LR and 
use units where h — ks — gfJ-B = e = 1. 

In order to calculate observable quantities for the sys- 
tem with Hamiltonian we find it convenient to use a 
fermionic representation of the local spin operator, 

& = 2 fj^Jl' f-y' i (3) 

77' 

with canonical fermion creation and annihilation opera- 
tors /t, / 7 , 7 =11, which allows a conventional diagram- 
matic perturbation theory in the coupling constant g. 
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Since the physical Hilbert space must have singly occu- 
pied states only, it is necessary to project out the empty 
and doubly occupied local states. This is done by in- 
troducing a chemical potential A regulating the charge 
Q = fyf-y- Picking out the contribution proportional 
to e~P x and taking the limit A — > oo, the constraint 
Q = 1 can be enforced (for a more detailed description 
of this method see I). 

We will use the Keldysh Green function method 
for nonequilibrium systems, following the notation of 
Ref. 0. Keldysh matrix propagators are defined as 

fl-(T£) w 

where G R ' A and G K are the retarded, advanced and 
Keldysh component Green functions, respectively. Spec- 
tral functions are found as A = i(G R — G A ), and the 
greater and lesser functions as 

q>/< = (QK ± qR t G Ay 2 ^ (5) 

The local conduction electron (ce) Green functions at 
the dot in the left and right leads, and the pseudo fermion 
(pf) Green function are denoted by G a ^ a and Q^ d , respec- 
tively, with lead index a — L, R, spin indices a, 7, and 
Keldysh indices a, b, c, d. A corresponding notation will 
be used for the self-energy S, and its imaginary part, the 
self-energy broadening, is denoted by T 7 = i(E R — T, A ). 
The interaction vertex has the following tensor structure 
in Keldysh space 

Kt = \(5a b Tl d + Tl b 5 cd ), (6) 

where a, b and c, d refer to pf, and ce-lines, respectively. 

Since we consider only nonequilibrium situations in a 
steady state, time translation invariance holds, and the 
single-particle Green functions depend only on one fre- 
quency. The bare pf spectral function is given by 

Ay(uj) = 27n5(cj + 7 5/2), (7) 

and the Keldysh component Green function is given as 

0f(w)=»A r (w)[2n 7 A(w)-l], (8) 
where n-y\(u}) denotes the pf distribution function, given 
byn-,x(w) = l/(e^ +x ^ T + 1) in thermal equilibrium. We 
shall also use the shorthand notation 

M 7A = 2n lX {uj) - 1. (9) 

Assuming a constant conduction electron density of 
states N(0) = 1/2D and a bandwidth 2D, the local ce 
spectral function takes the form 

A{u) = 2nN(0)8{D - \u\) (10) 

with the step function 9{x). The Keldysh component 
Green function in lead a is then given by 

G K{u)=-iA{u)t & ^(^^}, (11) 

assuming the electrons in each lead to be in thermal equi- 
librium. 



II. SPIN LEVEL BROADENING AND SPIN 
RELAXATION RATES 

The coupling of the local spin to the leads introduces a 
broadening of the Zeeman levels, which depends on tem- 
perature, magnetic field and bias voltage. In the pseudo 
fermion representation for the local spin, the broadening 
is given by the imaginary part of the pseudo fermion self- 
energy. This level broadening enters into the relaxation 
rates of both the transverse spin components (S X ,S V ), 
where it accounts for the loss of phase coherence, and the 
longitudinal spin component (S z ), where it describes the 
relaxation of the local magnetization following a change 
in the magnetic field. The observable spin relaxation 
rates, 1/T 2 and are defined through the broadening 

of the resonance poles in the transverse, and longitudi- 
nal dynamical spin susceptibilities, and their calculation 
requires vertex corrections to be included in a consistent 
way. 

Following a brief discussion of the pf self-energy broad- 
ening, we determine the renormalized ce-pf interaction 
vertex in a steady-state nonequilibrium situation. The 
resulting vertex functions are used to calculate the trans- 
verse dynamical spin susceptibility, and later, in Sec. lIIII 
they will serve as building blocks for a calculation of the 
conduction electron T-matrix. 



A. Pseudo Fermion Decay Rates 

In paper I (Ref. El), we determined the on-shell imag- 
inary part of the pseudo fermion self-energy, including 
leading logarithmic corrections. For the purpose of this 
paper, we will only need the second order rates, disre- 
garding logarithmic corrections. For T = one finds for 
< V < B: 

r T = \qIrV, (12) 

T l = T T + 27rg 2 B, (13) 
with Ag 2 = g\ L + g 2 RR + 2g\ R , whereas for V > B > 0: 

T T = ~gl R (3V-2B), (14) 
T x = T^+2ng 2 B. (15) 

Notice that in the presence of a finite magnetic field, only 
the upper spin-level, here corresponding to spin down, is 
broadened when V — 0, as one would expect from sim- 
ple phase-space considerations. Broadening of the lower 
spin-level (spin up) is due to virtual transitions to the 
upper spin-level and occurs only in higher orders in g. 

For comparison, we list also the thermal decay rate for 
V = B = 0: 

r u = 3nTg 2 . (16) 
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X = X • X 




FIG. 1: Vertex equation for the pf — ce vertex, including 
only non-logarithmic irreducible parts. 



B. Vertex Corrections 

Early worls22i2ii22i on the dynamical magnetic suscepti- 
bility of a single spin 1/2, demonstrated how self-energy, 
and vertex corrections combine to yield the transverse, 
and longitudinal relaxation rates I/T2 and 1/71. I n 
Ref. I20L the vertex corrections were determined in the 
approximation where the imaginary part of the pf self- 
energy, r, is much smaller than temperature. A simi- 
lar approach is possible out of equilibrium, where it is 
the finite voltage, rather than temperature, which deter- 
mines the abundance of (inter-lead) conduction electron 
particle-hole excitations. In the approximation where 
r « the dominant corrections to the ce-pf interac- 
tion vertex (Keldysh) tensor simplify substantially and 
the corresponding vertex matrix equation can be solved 
analytically. Since we assume that V 3> Tjc, perturba- 
tion theory is valid and r ~ g 2 V <C V is indeed a sound 
approximation. We shall consider the case where T -C V, 
which will best reveal the salient nonequilibrium features 
of the problem. 

To calculate vertex corrections and their interplay with 
self-energy diagrams, we have to solve the vertex equa- 
tion depicted in Fig. ^ We remind the reader that phys- 
ical quantities are proportional to e~^ x within our pro- 
jection scheme. Therefore we have to keep track of two 
contributions to the vertex 



A~ cd \ cd 1 A a cd 
ab = A a (, + A a b, 



(17) 



where °A is independent of A, and A A vanishes as e~@ x 
in the limit of A — * 00. We shall first determine °A and 
then, in a second step, A A. 



Voltage induced Particle-Hole Excitations 



We start by discussing the properties of the ce polariza- 
tion bubble in Keldysh space, which enters as one build- 
ing block in the vertex equation shown in Fig. ^ The 
convolution of two conduction electron Green functions 
has the greater component 



"• n <5>) = /^G>( £ + fi)G<( £ ), 



(18) 



and in general, the convolution of different Keldysh- 
components gives rise to the polarization tensor 



(19) 



It is convenient to form the contraction of this tensor 
with the exchange constants J aa i /4 at each end, and thus 
define an effective second order interaction by 
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nl c (n) = -j 2 a Xn d dc'm. 



(20) 



Contracting again this ce polarization tensor with two 
bare Keldysh vertices yields the pf interaction tensor 



D a'b a c! dTrd'c' \ cd' 

Db'a — A a 'ft lldc A a fe' 



(21) 



where the Langreth rules (cf. Ref. 23) have been em- 
ployed to work out the contractions 



mi = 2u K , 



U c / C = 2IT 4 , 



Uf c = 2U R , 



(22) 



using the notation 1 = 2 and 2 = 1 for the Keldysh 
indices. As for the single particle Green functions, we 
organize these components in a triangular matrix 



- - v 

and for <C D one finds that 



ii K {n) = ^ g 2 aa , {n + ^ a - (i a > ) coth 



(23) 



n 



4D In 2 
1 

7T 

Ma Ma 

~~2T 



(24) 



Notice that inter- lead particle-hole excitations do not sat- 
isfy the fluctuation-dissipation theorem as 



a a Il K {ty = coth 



~2T 



The lead-contracted polarization satisfies the following 
symmetries: 



n>/*/«M(_ft) = n<i K i A i R {n), 



(25) 



and for later use we quote the explicit formula for the 
greater component, U> = (U K + U R - II A )/2: 

n>(«) = ~ {gl R m + V)(l + N(Q + V)) (26) 

+(n - v){i + N(n - v))} + 2 ffo n(i + n(si))} , 

where N(fl) denotes the Bose-function and II > (il) = 
for il > 2D + V. In terms of this function, the second 
order pf decay rate may be written as 



T 7 (w) = 26 iY U > (-lj - -y'B/2), 
with 9 11 ' = 5 11 i + 2t 7 7 " 



(27) 
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2. Basic Approximations 



which is found to have the following components: 



The following calculations are based on self-consistent 
perturbation theory to order g 2 for self-energies and 
vertex corrections. As explained in detail in I, self- 
consistency is essential to obtain for example the correct 
magnetization out of equilibrium. For this paper, we 
need a self-consistent resummation of diagrams to inves- 
tigate how divergences of bare perturbation theory are 
cut off to lowest order in the interactions. 

However, for non-singular quantities like the lowest- 
order self-energy, self-consistency only gives rise to sub- 
leading corrections which we need not keep track of. For 
example, it is sufficient to approximate the retarded pf 
propagators (double-dashed lines in the diagrams) by 



■ 7 B/2 + ir 7 /2' 



(28) 



where T 7 , given in Eqs. (|12I15(I . denotes the on-shell de- 
cay rate calculated in bare perturbation theory. We ne- 
glect contributions from Re£ 7 (w) which can be absorbed 
in a redefinition of B and g, and which give rise only to 
subleading corrections in the following. 

To show, formally, that self-consistency does not 
change this result, one can use the fact that typical in- 
tegrals are dominated by integrations over frequencies in 
a window of width T around the Zeeman levels. Since 
the various Keldysh components of II vary slowly with 
frequency, i.e. [H(uj + T) — II(cj)]/n(ct;) ~ T/V, we may 
therefore use T/V ~ g 2 as a small expansion parameter 
(this will be shown explicitly for the vertex corrections 
below) . 



3. Summing up the Ladder 



To leading order in T/V ~ g 2 , the renormalized vertex 
satisfies the diagrammatic equation depicted in Fig. ^ 
This equation generates a series of ladder diagrams with 
dressed p/-legs and bare ce particle-hole propagators as 
rungs. Contributions from diagrams with crossed rungs 
we omit as being of order T/V. In appendix lAl we es- 
tablish this relative smallness explicitly for the crossed 
4'th order vertex correction and we expect higher order 
corrections to work in the same way. 

An iterative solution for the renormalized vertex starts 
with the attachment of two p/-propagators to the bare 
Keldysh vertex. This defines the tensor 



= 2A c J b ,G_TG 



b ar-fba 
7 —7' ' 



(29) 



"fired 
■y'V 12 



~f>V 21 



7'V 22 



X r R r K , 1 r R r R 

o c d y~/ yy + t ci1 y 1 t/y , 

0~cd {G^Gy + GjGy } 

^ T ld {G^Gy + G^Gy} , 
Scd GjGy + r\ d G^G. 



A 

7 "l'' 



(30) 
(31) 

(32) 
(33) 



One proceeds by attaching rungs, using the interaction 
tensor Bff£, and legs consisting of pairs of dressed pf 
propagators. This attachment consists of a contraction 
of Keldysh, and spin indices, together with an integration 
over the frequency circulating the individual sections of 
the ladder. To leading order in T/V, we may perform 
these integrals by neglecting the slow frequency depen- 
dence of the ce polarization functions compared to the 
rapid variations in the pf Green functions. Making use 
of the identity 



I 1 

a b 



1 



(34) 



products of Green functions may be expressed as either 



ft?(fi + w)#,(w) 



1 



n + (7-y)B/2 + i(r 7 + r 7 o/2 
i i 



lu + -/'B/2 - iTy/2 fi + cj + 7B/2 + ir 7 /2 



or 



G*(n + Lu)G5(uj) 



i 



n + (7 - Y)b/2 + i(r 7 - r 7 /)/2 
i 



u + j'B/2 + iTy/2 f7 + w + 7B/2 + T\/2 



and likewise for AR and AA products. Considered as 
an integral-kernel to be integrated with the various com- 
ponents of the polarization function, we may neglect the 
broadening and replace f2 by (7'— *y)B/2 inside the paren- 
theses in such products, and altogether this justifies the 
approximations 



&?(fi + w)£#(u;) 



2Tti5(uj + iB/2) 



q + (7 - i)B/2 + *(r 7 + ry )/2 ' 

0, (35) 



for a set of legs in the ladder. Notice that WalkesSfi has 
employed a similar approximation in the case of thermal 
equilibrium, utilizing the slow frequency dependence of 
the thermal ce-polarization. In this case, the RR and AA 
terms are neglected to leading order in T/T instead. 

Since the legs contain not only retarded and advanced, 
but also Keldysh-component Green functions, some of 
these loop-integrals will also involve the nonequilibrium 
^/-distribution functions n\(u>). This function is found 
by solving a quantum Boltzmann equation, obtained as 
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the Keldysh-component of the pf Dyson-equation with 
second order pf self-energies. Using the results of I, the 
solution at B — is found to be 

n\{uu) = n A (0)n < (w)/II > (w), (36) 

which, in the case where qlr ^ and T = 0, takes the 
form 



9i R (y-") 



n x (w) = n x (0) 



(9l L +g 2 RR )u+gl R (v+w) 



< LU < V 



9 2 lr(V-")-(9 2 ll+9rr)" ~V<UJ<0. 



(37) 

For T — > 0, n x (u) vanishes as V ^ T for u > V, 

and diverges as e~( w+v " T for w < —V. For < 
n x {u>) crudely resembles a Boltzmann distribution with 
T replaced by V/4. The distribution function clearly 
inherits the slow frequency dependence from II > and, 
to leading order in T/V, n x may therefore be treated 
as a constant, when integrated with the rapidly varying 
retarded and advanced pf Green functions. In the case 
of B > 0, the distribution function acquires a spin-index 
and the solution is generally more complicated (cf. I). 
However, the frequency dependence is still determined 
by II > , evaluated at arguments shifted by ±B/2, and 
therefore remains negligible. In either case, we are thus 
allowed to neglect the frequency dependence of n x , which 
renders Q K proportional to Q R — Q A by a constant and 
reduces all loop- integrals in the ladder to involve only the 
products (|35() or their complex conjugates. 

Omiting all RR and AA terms, V^ b d now simplifies to 



7 

4'V 12 



(38) 

i'v 12 = &cdG"G$, (39) 

W% = s cd {(i- M lX M YX )g*g*, - M lX M YX g*g$} 

+ri(M 7A -My A )^- (40) 
= -5 cd M lX gng$, (41) 

and performing the projection A — > oo, all p/-distribution 
functions vanish, i.e. M lX — > — 1, and we are left with 

#vgg(n + lo,lo) = -s cd T* b gX(n + w)g$(u). (42) 

Having performed the projection, it is now a simple 
matter to sum up the ladder solving the vertex equation. 
To keep matters simple we assume that B = 0, but once 
this special case is worked out, a generalization to B > 
will be straightforward. We begin by attaching the V- 
tensor (|42|l to the ce-polarization bubble defined in (|21|l . 
Working out the contraction, one finds that 



V c a d v {n + u J ,u J )B b % D (uj' -uj) = 



a'bi, I 



(43) 

-s cd TLG R (n + Lu)g A (uj)u> (<y - w ). (44) 

We should also attach the Pauli-matrices corresponding 
to the exchange vertices at the endpoint vertex and at 



each end of the polarization bubble. In zero magnetic 
field this yields the contraction 



7"'7" 7"7 ' j'f'" 'era' ' a' a 



-2r k , 



(45) 



which shows that the endpoint pf Pauli-matrix r k is car- 
ried through to the new external spin-indices. In this way, 
the Pauli-matrix at the endpoint vertex may be left out 
and the Keldysh vertex merely receives a factor of —2 
per rung. 

To second order in g, the vertex renormalizes to 

°A c a d b (n + u;\u;') (46) 
= Aft - f^VlUn + w, u)B$(u' - u) 



= 2 \ T cd S ab 



>cd 



'ab 



+ IT, 



„ 2n>(u>') 



n + ir 



where the left superscript is to remind us that the limit 
of A — * oo has been taken. The integral over to is per- 
formed using the ^-function from the iL4-product of pf 
Green functions and T is the spin-independent (B = 0) 
single pf self-energy broadening. 

Attaching a set of pf Green functions to this second 
order vertex correction, we notice that, after projection 
and discarding again all RR and AA products, we have 



a'b' 



c 3 nb' a nba! 



-V C{ 

v a 



cd 



(47) 



which in turn implies the fourth order correction 



0X C d( 4 ) 

= 2 



dej' iU>{uj') 



1 



2tt fl + iT 

2n>(o) 2n>(w") 



Vf$(Cl + tj',u/)B9*(u"-u/) 



2 5cd ^ aa n 



(48) 



From these two lowest order corrections it is clear how 
the further attachment to the ladder will generate a ge- 
ometric series, and the vertex function 



3 2ffl>(w) 



(49) 

therefore solves the diagrammatic equation in Fig^ in 
the limit A — > oo. We employ the suggestive shorthand 



r s = i(r t + r;) + r„ 

= 7Tg 2 LR V for B,T« V, 



(50) 



and, as will be demonstrated in the next section, this 
is indeed the spin-relaxation rate. In the present case 
of zero magnetic field and isotropic exchange couplings, 
the longitudinal, and transverse rates are identical and 
thus r s = T2 = T\. In the case of anisotropic ex- 
change couplings (or in the presence of a finite magnetic 
field), spin- flip, and non-spin- flip vertices receive different 
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corrections and the two rates become discernible. The 
anisotropic case will be discussed in Sec. I1VI The first 
term in (jf>0p arises from the self-energy, Eqs. (|12H15|) . 
the second one, T v — 2II > (0), is the vertex correction. 
Notice that only vertices with identical ingoing and out- 
going ce Keldysh indices are renormalized. 

So far, we have only determined the A — > oo limit of 
the vertex, but we need also the second contribution, A A 
in Eq. I|17fl . which is proportional to e _/3A . Having solved 
for °A already, we are left with the vertex equation 



x K d b (n + tS, w') = -2 / pA c J b „ (Cl + w,u) (51) 

J 27T 

xg b " a \n + u )g b ' a '\u)B&(u/ - u), 



which we find to be solved by 



4T, 



U_ ab (w) 



For c ^ d one obtains 



n 2 + r s 2 

[n>H(r Q 3 b + ir 2 ,) + n<(^)(r a 3 h - iT 2 ab )] 



(52) 



i[n x {Sl)-n x {G)]„a* 



rr», 



(53) 



which is neglected due to the slow frequency dependence 
of nx (ft). It is worth noting, however, that for B ^ this 
term will in fact be proportional to the magnetization 
and thus provide an important renormalization of the 
T cd t erm of the interaction tensor. 

This completes the solution of the vertex equation and 
we may now proceed to determine its influence on physi- 
cal observables. In doing so, one has to attach a pair of pf 
Green functions to the renormalized vertex, and most of- 
ten one may therefore continue to use the approximation 
(|35|l . Since the dependence of the vertex on the relative 
frequency ui is set by U_ ab (uj), one can safely set ui to 
and consider the vertex as a function of alone. With 
I\, = 2II > (0), the renormalized vertex then simplifies to 

A c d(n) = \rl d 5 ab + )-5 cd L ab (n), (54) 



where L ab =° L ab + x L ab , with 



iT v 




Sl+i(r+2I\.) 
Q+iT s 



n+<r f 



and 



x L ab {n) = -2n A (0) 



4r s n ab (o) 

n 2 + r? + '' ah n + ir s 



ah 

(55) 



(56) 



where n\(0) cx e ^ x . Using this result we can now calcu- 
late physical quantities like susceptibility and T-matrix. 



> 



FIG. 2: Dynamical susceptibility. Triangles refer to external 
measurement vertices. The black (emission) vertex is renor- 
malized like the interaction vertex in Fig. except that the 
two external ce legs are removed. The other (absorption) ver- 
tex remains undressed. 



C. Dynamical spin susceptibility 

In order to uncover the physical meaning of the rate T s 
introduced in Eq. I|5U|I , we include here a brief discussion 
of the transverse spin susceptibility: 



x x R (t)=m)([s-(t),s + (o)}). 



(57) 



The transverse spin relaxation rate, T%, is defined as the 
broadening of the resonance pole in this response func- 
tion, and as will be shown below, T s plays exactly this 
role. Throughout this Section, we may therefore use 
r s = T2- With a suitable generalization of r„, entering 
Eq. H50fl . which will be given in Sec. II VI this identifica- 
tion holds also for anisotropic coupling. 

Translating to the pseudo fermion representation on 
the Keldysh contour, the transverse susceptibility is cal- 
culated from 

^ X (r) = - i (- i ) 2 <T CK {/ i t (r)/ T (r)/ T t (0)/ i (0)}), (58) 

which in turn leads to the Feynman diagram in Fig. [21 
when including vertex, and pf self-energy corrections. 
The bare absorption, and emission vertices are given as 
lib = Tf^fc and lah = -j2 T ab> respectively (cf. I). The 
absorption vertex is kept undressed and the emission ver- 
tex renormalizes like the interaction vertex-component 
\/2A a l, whereby 



- x R (ty = i J ^Ail(n+cj,ij)g^(n+cj)gf(u;). (59) 



Notice that the canonical ensemble average, enforcing 
single occupancy on the dot, is carried out by dividing 
the A-dependent grand-canonical average by (Q}\ and 
taking the limit A — > oo (cf. I). This procedure affects 
only the pf distribution functions and allows to neglect 
all terms proportional to squares, or higher powers of 
n 7 \. Working out the contractions, we arrive at 
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2?r I 



>All(n + uj,w)2 [ ni (fj) - n^n + u)] gf{n + w)gf (w 



A A}} (ft + w) + A AM (ft + uj, uj) - 2n t (ui)°Ul (ft + uj,cu) Qf (ft + w)Sf (w) 

A A& (ft + w, w) - A AM (ft + w, w) + 2n T (ft + w) A|| (ft + w, w)] 0f (ft + w)SjV) } . (60) 



The important fact that the final result is proportional 
to n 7 is ensured by the relations 



°A n 



21 



-°A 1 1 { = A 2 1 1 1 - A^ = 1, 



(61) 



as such a constant drops after integrating over Q R Q R or 
G A G A . 

In the limit of S — ► 0, the factor of ni(uj) — n.f(ft + uS) 
in the first term is of order ft/V, and therefore we are 
forced to keep also the other terms involving Q R Q R or 
gAgA j n .j-j^g cas6j we nave to keep the full dependence 

of the vertex on two frequencies, but since for example 
the parts of the vertex which are retarded with respect to 
uj integrate to zero with R R , matters simplify substan- 
tially. The first square bracket can simply be replaced by 



2n x (0)n<(uj) - 2n> (uj)n x (uj) 



in\(u>) 



(62) 



and inserting now the nonequilibrium distribution func- 
tion given by (|36[) . the first two terms of this expression 
are seen to cancel. We emphasize the fact that this im- 
portant cancellation takes place only when using the cor- 
rect distribution function, i.e. the solution to the quan- 
tum Boltzmann equation corresponding to second order 
pf self-energies. 

The term involving Q A Q A works in a similar way, 
and using the approximation Q^"' A (Q + u>)Q^ , ' A (u>) ps 
— <9 w (w ± iO-f)" 1 , valid to leading order in max(|ft|,r)/V 
when integrated with the slowly varying distribution 
function, the last two terms in 16U|) may be evaluated by 
partial integration. The first term comes with a factor 

of + uj,Lu)g R (n + uj)g a (uj) « 2mS(w)/(n + %r 2 ), 

and altogether one finds that 



X 



M iT 2 



(63) 



for max(|0|,r) <§; V. The prefactor is independent of B 
and is obtained as the derivative — n'(0), with n(uj) given 
by (|37l) and with the replacement n\(0) — > 1/2, due to 
the normalization by (Q}\ before projection. The zero- 
frequency limit obeys X^O) = M/B, like in equilib- 
rium, and the nonequilibrium magnetization was found 
in I to be 



M 



(91l+9 2 rr + ^Ir)B 



(64) 



similar to a Curie-law with 1/T replaced by 4/V. Notice 
that the result llO.'ill . has been obtained also in Ref. llH 
using a Majorana-fermion representation. 

In the case of a finite magnetic field, the factor of 
nj,(w) — n-|-(Q + uj) in the first term of i|6U|) will be of 
order B/V. For B > max(|ft + B|,r), this term will 
therefore dominate the other terms involving Q R Q R or 
g A gA p Qr ^ > o, the vertex renormalization is modi- 
fied, but since we only need to consider the first term in 
(|60|l . only a single component is needed. For this partic- 
ular component the generalization is straightforward and 
one finds that 



°A&(fi + w,u;) 



2iTi>(uj- B/2) 
Q + B + iT 2 



(65) 



where T 2 is given in Eq. I|5U|I and depends now on both 
V and B (see Eq. (I68|) below). The integral over uj is 
performed using the approximation l|35|) , and the suscep- 
tibility is found to be 



~X 



«(n) 



M 



n + B + iTo 



(66) 



valid for max(|ft + B\,T) < min(B, V) 

In the intermediate regime where B <C min(|ft+.B|, T), 
one would need to generalize also the A-dependent part 
of the vertex to the case of B > 0. However, we expect 
that cancellations, similar to those found in terms like 
(|62|l at zero field, will take place also at finite B, once 
the correct B-dependent distribution function is used. 
In this manner, we expect the general formula for the 
susceptibility to be simply 



M B + iT 2 
~B ft + B + iT 2 



(67) 



valid for max(|ft + B\,T) <C V. This function obviously 
has the correct asymptotic behaviors, corresponding to 
(|63[) and (|66ll . and is consistent with the equilibrium 
resullM. 

For completeness, we state here the relevant asymp- 
totics of r 2 as a function of V, B and T: 



max(T, B) < V 



r-- { <9 2 ll + 9 2 rr)B/± , max(T, V) « B 

<9ll + 9rr + ^9 2 lr)T , max(B, V) « T 



(68) 
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In the equilibrium limit, V = 0, this corresponds 
to the result obtained in Refs. |53 and l2ll T 2 ~ 
7T.g 2 max(T, B/A). 

III. CONDUCTION ELECTRON T-MATRIX 

With the renormalized vertex at hand, we now proceed 
to calculate the conduction electron T-matrix, including 
the leading logarithmic corrections. The T-matrix, T aa i, 
is of great physical significance, insofar as it describes the 
scattering of conduction electrons from lead a' to lead 
a, and thereby also the transport across the dot. It is 
determined from the conduction electron Green function: 

(69) 

In cases where the exchange-tunneling Hamiltonian is 
derived from an underlying Anderson model, i.e. from a 
single quantum dot in the Coulomb-blockade regime (cf. 
e.g. Ref. ll3jk one has J\ R — JllJrr and only one of the 
eigenvalues of the 2x2 matrix T aa i is finite. In such a 

I 

Im[T* Q (fi)] = 
Im[T£ R (fi)] - 

with = — /in = V/2. Within bare perturbation the- 
ory, the T-matrix diverges close to each Fermi surface, or 
more precisely, for f7 — > Ma, some of the logarithms are 
cut off by the voltage V = [i^ — /i^ while others remains 
unaffected. In this sense voltage and temperature act 
very differently as T would cut off all logarithmic terms 
uniformly. The central question formulated in the intro- 
duction is, how the logarithmic divergences which remain 
for T — > and large V are cut off when the perturbation 
theory is properly resummed. To find the correct cut-off 
to order g 2 , we have to replace the bare Green functions 
and bare vertices in Fig. E]by the dressed ones. 

As the second-order diagram in Fig. Ogives only a finite 
contribution Im[Tf Q ,(fi)] = - J2 a " 9aa"9a" a ', the 
inclusion of self-energy, and vertex corrections will pro- 

I 




FIG. 3: Diagrams for the conduction electron T-matrix, with 
dressed pf propagators and dressed interaction vertices (black 
dots). 



situation, Im[T QQ / CT (u;)] is, at low energies, directly pro- 
portional to the s pec tral function of the electrons on the 
dot (see e.g. Ref. yjand references therein). This spec- 
tral function can be measured directly by tunneling into 
the dotS, and henceforth we shall focus on the imaginary 
part of T aa : . 

In Fig. |3| we show the two diagrams contributing to 
the T-matrix to third order. Within bare perturbation 
theory (i.e. using bare vertices and Green functions in 
Fig. |3J), one obtains the following intra- and inter- lead 
components at T, B = 0: 

(70) 
(71) 



duce only subleading corrections of order g , as can be 
shown by an explicit calculation. 

The fate of the logarithms arising to order g 3 is more 
interesting, and in the following we will therefore care- 
fully evaluate the second diagram in Fig. [3] This contri- 
bution involves the spin-contractions 

T ~f7' T y'j" T j"~f T a<7" T a"a' T cr'cr = (^) 

for the Peierls, and 

t^t^t^t^^tX = -24 (73) 

for the Cooper-channel. Writing out the sum of these two 
types of diagrams, corresponding to different orientations 
of the pf -loop, one finds that 



3tt 
16iV(0) 

16iV(0) 



(sL + 9lr) 



1 + 2g aa In 



D 



\Q - fi a 



9lr \ 2g, 



1 + 2g LL In 



D 



|0 - fi a \ 



r- ^9d9 2 LR In 
2(<7fl fl + 9lr) ln 



D 



D 



|«-Mfl| 



r 



x {A$ (w, w + e)g ba MAfi (cj + s', u)g_ b ' a ' {w + e ')A^ (w + e, u> + e')G b " a " (u + e) 
-A^( W , W - e')S b »A^,> - e,u;)g b " a " ( w - e)A^',(c - e',u - e)£ hV (c - e ')} ■ (74) 
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This Keldysh contraction has a total of 256 terms, of 
which only a few will contribute in the end. Some will 
involve a C/ 21 , which is zero, and others will involve a 
product of more than one lesser-component Q < 1 which, 
being proportional to higher powers of the pf distribution 
function, will vanish faster than (Q)\. Since the Keldysh 
representation contains Q < as part of Q K = Q K + Q y , 
it is a daunting task to isolate all contractions with only 
one factor of Q K . Nevertheless, since we are dealing here 
with a trace over the pf Keldysh indices, we are free to 
work in a more convenient basis for the pseudo fermions. 
Thus choosing 



where G^(e) = l/(s-iT s ) are Green functions broadened 
by T s rather than T/2 and we use the shorthand notation 
Jaa"a'"a' = J aa" J a" a'" J a '" a ' ■ Already at this stage, it 
is apparent that the vertex corrections have served to re- 
place twice the pf self-energy broadening by T s . Making 
use of the basic integrals, 



D 2 



D sgnje + a)(e + b) 

D (e + 6) 2 + r2 V(&-a) 2 + rf 



(81) 



and 



c 



1 1 
o 1 



i -i 

1 



(75) 



the Keldysh matrix Green functions may be transformed 
as 



Q-CGC =( gA 



(76) 



which has the nice property that Q becomes diagonal 
after projection. The renormalized vertices may be con- 
sidered as functions of only one frequency and therefore 
take the form 1)54)1 . which we write loosely as A = S + L. 
For opposite ce Keldysh-indices the vertex retains the 
structure of the identity-matrix S under the transforma- 
tion. For equal ce-indices, the matrix L a ^ transforms to 



L = (LC 1 = °L + X L, 



(77) 



D 



d (e + b) 2 + r s 2 



(82) 



representing a broadened logarithm and a broadened 
sign-function, respectively, the remaining integrals over 
e and e' are straightforward. 

The first line of the integral (|8U[) involves a convolution 
of G K with Q A , which yields 



-iN(0) < In 



D 2 



(ci- f i a „) 2 + r 2 



2i tan 



x /n - Ho 



D 



This term is multiplied by the convolution of G^> with 
Qy, equal to iG^:,(fl + iT), and altogether the first line 
yields the imaginary part 



where 



°L 



1 

<P -1 



(78) 



-2n A |j^^AT(0) 2 ln( (o _ M ^ 2 + r , 



(83) 



and 



2n A (0) 



£1 

2 



4r s n"(o) 



2ir„ 



Vil 

2 



4r,n"(o) 
o 2 +r 2 



(79) 



with ^ = ^ and ^ = n^r B - 

In this representation the contraction in Eq. 174f) be- 
comes manageable and one has to deal with merely 8 
different types of terms. The full contraction is worked 
out in Appendix iBl resulting in 



(80) 



de f de' 

lb J Zir J In 

-[G&,(n + E )G$„(n + e / ) 

+G%, (fi + e)G£,„ (fi + e ')] Qi {e)Gl (e ~ e') } , 



Using a spectral representation for the ce Green func- 
tions, the remaining two lines of l|80|l can be brought to 
the form 



2n x 



sgn(e - fi a i,i) 



(u> — — iT s )(u> — e — iT s ) 

sgn(e - (j, a „) 

' (e ~ n - iT s )(uj - e + iT s ) 



(84) 



The u> integral in first term vanishes in the limit D — * oo, 
and keeping D finite this term remains smaller than the 
second term by a factor of Q./D or T s /D. Keeping only 
the second term, the imaginary part takes exactly the 
same form as (|83H . and finally we obtain after projection 
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Im [T« (fl)] 



16JV(0) ^ 



1 + g aa In 



_D 2 



(n- Ma )2 + r 2 
i 



2g 2 LR9d In 



D 2 



(n + MQ )2 + r2, 



(85) 



with no summation over a implied. We find precisely 
the result of bare perturbation theory, Eq. H7U|) , but now 
with the logarithmic divergences cut off by T s . The same 
conclusion also holds for Tlb- This is the central result 
of this paper. 



IV. ANISOTROPIC COUPLINGS: Ti VS. T 2 

The longitudinal, and the transverse spin-relaxation 
rates, 1/Ti and I/T2, have rather different physical in- 
terpretations. It is therefore interesting to determine 
which combination of the two rates actually controls the 
logarithmic divergences. In the previous chapter we re- 
stricted ourselves to the case of zero magnetic field and 
isotropic couplings, and in this case we cannot distinguish 
between the two rates as 1/T% = I/T2 = T s . 

To discriminate between the two rates, even for B = 
0, we generalize the exchange interaction to involve two 
different couplings 



7 7 o o 



T 2 , ■ 

7 7 



- 2 > ) 



Z J 
■' < v 



3 3 
1 I'l cr'cr ' 



(86) 



and we may now repeat all calculations above, keeping 
track of separate spin-flip and non-spin-flip processes. 
Since we consider only the case of zero magnetic field, 
the pf self-energy broadening remains spin-independent 
and we obtain from Eq. $2$, for V > T and B = 0, 



T = T^ =T l = -( z g 2 LR + 2^g 2 LR )V. 



(87) 



The vertex corrections now take a different form, de- 
pending on whether or not the spin is flipped at the 
vertex. For T = B = and finite V we obtain 
= n z gj jR V/4: for the spin-flip vertex and T z = 
ir{2 ± g 2 LR — z g 2 LR )V/'i in the case of no spin-flip. There- 
fore, the longitudinal, and the transverse spin-relaxation 
rates are given by 



± r =T l =T + Tl=-K^gl R V, 



and 



7f — r 2 
J2 



T + Tt 



(88) 



(89) 



Notice that 1/Ti = for L g = 0. This is due to a can- 
cellation of vertex, and self-energy corrections, reflecting 
the conservation of S z in this case. 

How do these spin-relaxation rates modify the loga- 
rithmic divergences? A close inspection of the Keldysh 
contractions and the integrals carried out in Appendix IBl 
reveals that only the °A-part of the renormalized vertex 
connecting to the out-going ce-line (i.e. the left most ver- 
tex in Fig. |3 gives rise to a logarithmic divergence, and 
furthermore determines whether this logarithm is cut off 
by T2 or Ti depending on whether this vertex involves 
a spin-flip or not. Therefore, in the case of anisotropic 
couplings, Eq. I|85|) generalizes to 



Im[Tl,(0)] 



16A(0) 



Qol"o.' Qo."a.'" Qo.'"ol' In 



D 



2 g a "a' + ( Z 9a"a'" 9a'"a' + 9a"a"' Z 9a" / a') m 
I 



D 2 



(90) 



(n-(i a „) 2 + r 2 



Roughly speaking, two thirds of the logarithms are 
broadened by T 2 and one third by Ti. 

How are these results modified beyond lowest order 
perturbation theory? In Appendix ICl we investigate this 
question in the limit ^-g for finite z g. In this limit, 
the logarithmic singularities in correlation functions like 
(S~ S + ) resum in equilibrium to power-laws with expo- 
nents depending on z g. In Appendix IU1 we use a mapping 
of our problem to a non-equilibrium X-ray edge problem 
together with results by Ng2& and others^Zi 2 ^ to investi- 



gate how these power-law singularities are affected by a 
finite bias voltage and the associated current. We find 
that all these power-laws are cut off by a rate related to 
I/T2. This has a simple interpretation: for finite Z J a 
finite current is flowing through the system and the cor- 
responding noise prohibits the coherence of the two ex- 
ternal spin-flips at low energy. Close inspection reveals, 
that the second logarithm in Eq. I|90|l is calculated from a 
correlation function of the type discussed in Appendix IH1 
The non-perturbative results of the Appendix therefore 
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confirm our perturbative Eq. (|90|l . The first logarithm in 
Eq. (|90l) . however, arises from a different correlator (as 
one external vertex involves S z ) which we have not tried 
to calculate to higher orders in z g. 

Also in the presence of a magnetic field the situation 
is more complex and at present we do not know which 
combination of relaxation rates controls the logarithmic 
divergences arising for V ~ B. The vertex corrections 
depend on B and, as mentioned in Sec. Ill B 31 also the 
r^-part of the vertex renormalizes in this case. Further- 
more, the non-spin-flip vertex depends on the orientation 
of the incoming spin, and its two different components arc 
found only after solving two coupled vertex equations (cf. 
e.g. Ref.Ej. 

In many physical situations, \/T\ and I/T2 differ only 
by a numerical prefactor of order 1 and such a factor in 
the argument of the logarithms is not important. In this 
situations it is not necessary to keep track of differences 
of 1/Ti and I/T2, if one is interested in a calculation to 
leading order in 1/ In [max (V, B)/T K ] (cf. e.g. Ref.0). 



in order g 2 , to the number of conduction electrons passing 
the constriction per unit time (the factor of proportion- 
ality depends, however, on details of the model, such as 
e.g. anisotropics of J). We therefore interpret the sub- 
sequent attenuation of the Kondo effect as decoherence 
due to current-induced noise. 

Most formulations of perturbative renormalization 
group in equilibrium completely neglect the role of de- 
coherence and noise and focus instead on the flow of cou- 
pling constants. This is justified, as the typical rates are 
often much smaller than temperature T, which serves as 
the relevant infrared cutoff. However, since this is not 
the case in a nonequilibrium situation, decoherence has 
to be an essential ingredient in any formulation of pertur- 
bative renormalization group valid out of equilibrium^. 
We hope that our perturbative calculation, demonstrat- 
ing how this happens in detail, can serve as a starting 
point for future developments in this direction. 

Acknowledgments 



V. DISCUSSION 

In this paper we have addressed the question how, far 
out of equilibrium, the presence of a sufficiently large cur- 
rent prohibits the coherent spin-flips necessary for the de- 
velopment of the Kondo effect. In an explicit calculation, 
we have confirmed the expected answe»Siii*iSiii that the 
spin-relaxation rate cuts off the logarithmic corrections 
of perturbation theory. This implies that for T a Tk 
(i.e. for V » Tk, see Refs.00), the Kondo-model stays 
in the perturbative regime, which allows calculating its 
properties in a controlled way using perturbative renor- 
malization groupS. 

We have worked out this scenario explicitly for the 
imaginary part of the conduction electron T-matrix, tak- 
ing into account the joint effect of self-energy, and ver- 
tex corrections. In the limit of zero temperature and 
\h(V/Tk) 1, perturbation theory remains valid and 
the vertex corrections were determined by summing up 
diagrams to leading order in T/V ~ g 2 . Within bare 
perturbation theory, the T-matrix exhibits logarithmic 
divergences at the Fermi energies of the left, and the 
right lead, and we have demonstrated explicitly that the 
joint effect of dressing pf Green functions as well as ex- 
change vertices with voltage induced particle-hole exci- 
tations works to cut off these logarithms by T s — ^g\^y ■ 
Under certain conditions, the T-matrix can be identified 
with the spectral function on the quantum dot, which 
can be measured directly by tunneling into the do^2^. 

To reveal the physical significance of this rate, we have 
calculated the dynamical transverse spin susceptibility 
in the presence of a finite bias-voltage. This served to 
demonstrate that T s is indeed the spin-relaxation rate, 
broadening the resonance pole at uj ~ B in this correla- 
tion function. T s arises from the stirring up of inter-lead 
particle-hole excitations, and is found to be proportional, 
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APPENDIX A: VERTEX CORRECTIONS FROM 
CROSSED RUNGS 

In this Appendix, we argue that crossed diagrams of 
the form shown in Fig. 0] lead only to subleading cor- 
rections. In such diagrams, for example, the simple fre- 
quency dependence observed in Sec. Ill B 31 is no longer 
valid. We therefore evaluate explicitly the crossed 4th 
order correction depicted in Fig 0] and compare it to 14811 . 

The Feynman rules give the same prefactors in this 
case, and the contraction of spins yields 

(TVrV n T n ) y7 Tr[T m T J ]Tr[TV 1 ] = 20 r*, 7 , (Al) 

as opposed to 4r*/ obtained in the ladder-type correc- 
tion. The Keldysh contraction may be expressed in terms 
of the previously defined tensors V and B as 

v c J* (o + w, w)bW\u? - uj)g b, " a "(n + «') 

x G b " a "\u J + to" - cj')B&,,(u" - J), (A2) 

and using the identity Ij42[) . this may be worked out to 
give 

x 2n>( w " - u')g R (n + u')g R {u' - J 1 ). (A3) 
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FIG. 4: Vertex correction from crossed particle-hole excita- 
tions. Such contributions are smaller than the ladder-type 
corrections by a factor of T/V and are therefore neglected. 



This result is reached after eliminating a number of 
contractions which vanish under projection, that is, all 
terms proportional to one or more factors of G < , and 
apart from an overall factor of 5, coming from the differ- 
ent spin sum, this contribution looks very much like the 
ladder- type correction Q48)). However, the second pair of 
pf Green functions have the structure of an RR product 
as a function of u>', and when integrated with n > this 
makes (|A3J| smaller than (|48f) by a factor T/V. 

Notice that, in contrast to the ladder diagrams, the 
crossed diagram in Fig. 0] involves a loop-integral over 
QQQTi, which does not warrant the omission of RR and 
AA terms leading to (|42|) . However, keeping all terms in 
the ^-tensor, a rather lengthy contraction leads to a re- 
sult which differs somewhat from l|A3p . but nevertheless 
remains smaller than l|48|l by a factor T/V. We expect 
higher order corrections to work in the same way, which 
renders the ladder-type corrections dominant in T/V to 
all orders in the coupling. 



APPENDIX B: CONTRACTIONS FOR T R(3) 

In this Appendix, we work out the contraction of 
Keldysh indices in Eq. (|74|l . There is a total of 9 dif- 
ferent non-zero contraction of ce Keldysh-indices, each 
of which involve renormalization of either zero, one, two 
or all three vertices. This gives rise to a total of 2 3 = 8 
different types of pf traces, which we need to work out. 
If the two ce Keldysh indices are different, a vertex con- 
tributes with a factor of 5 a b rather than L a \,. Thus a 
term with all three vertices renormalized contributes with 
Tt[LQLQLQ], whereas a term with no vertices renormal- 
ized contributes Tr[SGSGSG]- Our strategy will be to per- 
form the contraction and the loop- integral over uj without 
including the A-dependent part of the vertex. After this 
has been done, it will be a simple matter to include the 
additional effects of °L, by going through very similar 
steps once more. 

We begin by listing a few useful facts about the rele- 
vant matrix products: 

°^=(|?« 2<f>g<-g A )> 6C= \ G o %) (B1) 



and 



Tr 



ai 6i 
ci di 



a 2 b 2 
C2 d 2 



«3 b 3 
C3 d 3 



a 1 a 2 a 3 + bic 2 a 3 + a x b 2 c 3 + b\d 2 c 3 + 

C1O2&3 + dic 2 b 3 + C\b 2 d 3 + d\d 2 d 3 . (B2) 

The lesser component Green function takes the form 
g < = ri\(Q A — Q R ), and neglecting their slow frequency 
dependence we may consider the pf distribution func- 
tions as constant prefactors. This allows us to expand 
all terms in products of three Green functions which 
are either retarded or advanced, and to use rules like 
Q\Q%Qz = QfG 2 Gi = 0, implied by the subsequent 
loop integration which can now be performed by closing 
in the half-plane with no poles. Notice that including 
the frequency dependence in either factors of n\ or n ab , 
coming from either propagators or vertices, would render 
such loop-integrals non-zero. Nevertheless, these contri- 
butions will be smaller than the terms which we retain 
by a factor T/V and can therefore be neglected. Further- 
more, the projection allows us to neglect terms which are 
proportional to Q < Q < or Q < Q < Q < . 

With these few rules at hand one may work out the 
following catalog: 



Tr 



( Q ^)i( ^) 2 ( ^) J 



(B3) 



An x {^{g R g R g A ~g R g A g A ) 
+(f>2(g A g 2 g 3 - gts 2 s A ) 



(t>3(grg 2 gf - <?i g 2 g 3 )} , 



Tr 



Tr 



Tr 



CmiCWMsgh 



4m 



( ^)i(^) 2 (°^) 3 



4n> 



(^)i( ^) 2 (°^) 3 



hg?g£gt 

Ar.Rn.A 



4n A {<i> 2 g? g?g$ 



^g^g^gi) 
hgtg^gz} 



The remaining four possibilities all vanish, and we are left 
with contributions from terms with either two or three 
vertices renormalized. Working out the loop-integral over 
uj, we get e.g. 



— 01^1 y 2 y 3 



dui 
2^ 



<l>(-e)g R (w) 



xg R {oj + e')g A (co + e) 
= l g A 2 (e)g A (e-e<) 7 (B4) 

where we have introduced the notation g A o (e) = (e — 
iFs)" 1 , and g A (e) = (e — ir) _1 for the double-broadened 
pf Green functions. We see that the vertex corrections 
serve to replace T by T s in products of certain internal 
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Green functions, and working out all the integrals, we 
obtain the following list for the Pcicrls-channel: 



— <f>iG R G R G£ 

Z7T 



2^ 



>2^i y 2 ^3 



2^ 



— <j>zG R G A G R 



hG A G A G R 



iGf 2 (e)Gf(e-e'), 

-iG^ 2 (e)G^(s'), 

l G?(e)G? 2 (e'), 

iG? 2 (e')Gf(e-e'), 

-iG£(e')G? 2 (e-e'), 

-iG?(e)G^{e-e'). (B5) 



As may be seen from Eq. <|74[1 . the corresponding prod- 
ucts for the Cooper-channel can be obtained from these 
by the shift of variables e — ► — e', and e' — ► — e. Using 
the fact that G R {—s) = —G A {s), one readily obtains the 
following list, to be used for the Cooper-channel: 



J^GMGi 



du> 
2^ 



/ % hG A G 2 R G : 

J 2^ fa^^G- 



i <->Br>A<->R 

>3feq y 2 »3 



2^ 

— faG A G A G R 

Z7T 



-iGP 2 (s')G£(e-e>), 
-iG? 2 (s')G?(e), 

iG A (e')G^M 
-iG£ 2 (s)G£(e-e'), 

iG?(e)G? 2 (e-e'), 

iG£{e')G? 2 {e-e'). 



(B6) 



It is now straightforward to carry out the contraction 
of ce Keldysh indices in Ea. H74|) . and one finds the com- 
bination 



G R G R Tr[( LG)i( LGM°LGh 
+G K G R Tr[(°LG) 1 (°LGH6Gh 
+ (G R G K + G k G a )Tt [( q LG)i(6G)2(°LG)3 

for the Peierls, and 

G r G r Tt 



( LG)i( LG)2(°LG) 3 

( Q LG)i(°LGMSGh 
(G R G K + G k G a )Tt \(5G)i( LG) 2 ( LGh 



+G K G R Tr 



tribute the integral 

^ (3) (n) = ^ w /g/g (B7) 

x {G%,(Sl + s)G R „(n + s')G£(e)Sr(e') 
- [G R „(n + e)G%„(n + e') 

+Gg,(fl + e)G A „(Q + e')] G£(s)G£(s -a')}. 

To include the effects of A L, one may go through the 
same steps and build up a similar catalog of terms. We 
have to include all terms with exactly one factor of X L, 
since terms with two or three factors vanish faster than 
(Q)\ under projection. To leading order in T/V, there 
will still only be contributions with either two or three 
vertices renormalized. Whereas L ended up contributing 
only with its 21-entry, <f), this entry is zero in L and 
instead one finds only contributions from its 12-entry, tp. 
A typical contribution from the Peierls-channel now takes 
the form 



duo r 
2^ 



-Tr 



(°LGM x LGH6Gh 



= 2n x 



4n> 



doj 
to* 



GffaGiGi 



G£{e')Gi{e), 



for the Cooper-channel. Together, the two channels con- 



and a term like this eventually adds up with a sim- 
ilar term from Tr[(°LG)i{°LG)2(SG)3\, having a 1 in 
place of the factor of iT v /(e' — iT s ), to contribute 
in\G A B {e')G A s {£)- Working out the full contribution, 
from both the Peierls, and the Cooper-channel, one finds 
that all surviving terms combine in similar ways, and the 
total effect of including X L is therefore simply to replace 
T by T s in {B7|). This finally leads to the integral l|gOf) 
quoted in the main text. 



APPENDIX C: CUTTING OFF X-RAY EDGE 
SINGULARITIES IN THE ANISOTROPIC 
KONDO MODEL 

In this Appendix, we investigate the anisotropic Kondo 
model in the case of a vanishing spin-flip coupling J = 
and finite Z J . In this limit, certain equilibrium correlation 
functions are singular at the Fermi energy, they display 
the so-called X-ray edge singularities whenever the spin 
is flipped. In the following, we investigate how these 
singularities are modified in the case of a finite voltage. 

Even for J = a finite current is flowing through 
the system as z Jlr ^ and we therefore expect that the 
associated noise will cut off all singularities. Fortunately, 
a very similar pro blem has been solved exactly by Ng2£ 
(see also Refs. W\ and 12^1 . who considered the effects of 
suddenly switching on the tunneling between two (non- 
interacting) leads. 

We will show, that our problem (for J = 0) can be 
mapped exactly on the one solved by Ng. The fact that 
this is possible is not obvious as he considered a situation 
where for times t < ti no current is flowing, whereas in 
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our case the same current passes the dot before and after 
the spin-flip . 

Ng considered the Hamiltonian 26 

H X = H Q (V) + V r a / a c| rtt ,Ca k ff(t/-t)e(t-t i ) l 



,k,k' 



(CI) 

where H (V) = J2 a ,k,v( e * - Ma)c^ k(T c Qk(T describes the 
two leads with the bias voltage V = /i£ — [ir. The tun- 
neling between the left and the right lead (and a poten- 
tial scattering) is switched on for times between ti and 
t j . This generalization of the usual X-ray edge problem 
to two different Fermi seas was solved by Ng2&, using 
a generalization of the method devised by Nozieres and 
De Dominicis for the problem with only a single Fermi 
sea. He finds that the relevant spectral function exhibits 
power-law singularities near each of the two Fermi ener- 
gies in the left and right leads, which are, however, cut 
off by a voltage induced broadening given in terms of 
complex phase-shifts (see Ref. Efil for details) 5l/r.i by 



V 
2^ 



Im[5 L 



(C2) 



For -Lg = and B = 0, the Kondo Hamiltonian 
reduces to two separate potential scattering problems for 
conduction electrons of spin up and down, respectively, 

H = H (V)+ ( Z ^'^72)ci, kV ,r a 3 v c akCT . (C3) 

a. a' ,k,k' .a. a' 

and we want to study the effect of a single spin- 
flip, i.e. correlation functions like (S~(tf)S + (ti)) or 
([(4, T (i/)£ - (f/)] [S + (U)c a ^(U)}) (which is related to the 
T- matrix). For these correlation functions, the spin 
points down for t < ti, i.e. S z = —1/2 and H(t < ti) — 
H t = H (V) - \ E a , a ',k,k',^ z -4'aC Q , /k/0 .,T^ CT c O!ktT . To 



map Eq. (|C3j) onto Eq. I|C1|1 we note that S z 
ti < t < tf and therefore 



1/2 for 



H = H i+ £ ('Jctaffl^^CatoOitf-tWt-U). 
a,a',k.k',cr,(7' 

(C4) 

Hi can easily be diagonalized in terms of scattering 
states. Scattering states coming from the left (right) 
lead are occupied according to the left (right) chemical 
potential and therefore, (|C4J| takes the form (|C1|I when 
rewritten in terms of those scattering states. 

To determine the scattering states of Hi, we represent 
for convenience the two semi-infinite leads by infinite chi- 
ral wires of right-movers. In this representation, the scat- 
tering wave-functions ^^"(x) describe the amplitude of 



plane waves coming from lead a 



$&(x) = [6{-x)5 a , a + 6{x)S a>a ] e lkx , 



(C5) 



where x < (x > 0) refers to incoming (outgoing) waves 
in lead a' . The scattering matrix S a i a is determined from 
the Schrodingcr equation 



-iv F d x 8 a . a „ - - z J a , a »5(x)\ $^ Q (x) = e^(x), 

(C6) 

and regularizing the delta-function by using 8(0) = 1/2 
we obtain 



Sll — 



~>LR 



1 - ( z 9 2 lr ~ z 9ll z 9rr)/M + io[*g LL - z g RR )/8 
1 + ( z 9 2 lr ~ z 9ll z 9rr)/^ ~ i°( z 9LL + z 9rr)/^ 

z 9lr^/8 

1 + ( z 9 2 lr ~ z 9ll z 9rr)/M ~ M z 9ll + z 9rr)/& 

(C7) 

with z g aa , = N(0) z J aa , and N{0) = l/v F . 

Rewriting Eq. I|C4(I in terms of these scattering states, 
we can read off the potential in Eq. (|C1|) 



v a , a = -J2 Z J >0 [< a '(o)V ^"(0). (C8) 

w 

Using this formula and the results by Ng2&, one can easily 
work out the relevant correlation functions when taking 
into account that the spin-up and spin-down problems 
separate. The corresponding correlation functions are 
therefore multiplied in the time-domain and convoluted 
as a function of frequency. We will not display the rather 
lengthy formulas, but only note that all divergences close 
to the two Fermi levels are cut off by the appropriate 
relaxation rates ljC2fl [the rates for spin-up and spin-down 
add as e - r T*e- r i' = e-( r T+ r i)*]. 

To make contact with our perturbative results, we 
will now consider the case of small Z J . In this limit 
V a ' a « Vjjip/2. Inserting this into Eqs. (lid) and (llf) 
of Ref. |2y, determining the complex phase-shifts 8l/r, 
expanding the result to leading order in V a > a and adding 
spin-up and spin-down contributions, we find 



r x = 2 V \ z 9lr\ 2 , 



(C9) 



which coincides with our F2 = I/T2 in Eq. Q89)). in the 
limit of g — > 0. Note that the first logarithm in Eq. I|9U|) 
arises from a diagram with S z at an external vertex. 
Therefore the corresponding correlator is not of the X-ray 
edge form discussed in this Appendix. 
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